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We investigate polariton-polariton interactions in a semiconductor microcavity through two- 
dimensional Fourier transform (2DFT) spectroscopy. We observe, in addition to the lower-lower and 
the upper-upper polariton self-interaction, a lower-upper cross-interaction. This appears as sepa¬ 
rated peaks in the on-diagonal and off-diagonal part of 2DFT spectra. Moreover, we elucidate the 
role of the polariton dispersion through a hne structure in the 2DFT spectrum. Simulations, based 
on lower-upper polariton basis Gross-Pitaevskii equations including both self and cross-interactions, 
result in a 2DFT spectra in qualitative agreement with experiments. 
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I. INTRODUCTION 

The strong coupling between quantum well excitons and 
photons confined in a microcavity gives rise to two new 
eigenstates: lower and upper-polaritons. Furthermore, 
polariton-polariton interactions (anharmonicities), medi¬ 
ated by the nonlinear interaction of excitons, provide a 
wide range of rich physics. In fact, the lower and upper 
polariton states are no more exact eigenstates, because 
they are defined in a non-interacting (harmonic) exciton 
system [T]. The exciton-exciton interaction introduces 
not only lower and upper-polariton self-interactions but 
also lower-upper cross-interactions. In quantum chem¬ 
istry, this type of problem is known as “normal mode 
versus local mode” problem [2]-[4]. In our context, the lo¬ 
cal and normal modes are respectively the exciton-photon 
and lower-upper polariton states. Although a wide range 
of research has been made on lower polaritons, such as 
superfluity [5] and its Bogoliubov excitation spectrum 
mn\, the self and cross-interaction between lower and 
upper-polaritons has not yet been fully investigated. In 
fact, it is difficult to distinguish lower-upper polariton 
cross-interaction from self-interactions through conven¬ 
tional one-dimensional (ID) spectroscopy [ 8 ]. In order 
to enlighten the polariton interactions, it is useful to em¬ 
ploy a two-dimensional Fourier transform (2DFT) spec¬ 
troscopy technique. 

2DFT spectroscopy is a powerful tool to investigate 
coherent couplings and vibrational anharmonicities of 
molecular vibrational states ElEl- One advantage of a 
2D spectrum is that we can associate each peak of the 
spectrum with different Liouville-space pathway through 
double-sided Feynman diagrams mm- With this idea, 
we can identify dominant nonlinear optical pathways in 
four-wave mixing (FWM) signals, which we cannot access 
with conventional ID spectra. This method has been 
extended to investigations of electron-hole many-body 
properties in semiconductor systems mm- Those re¬ 


searches revealed the importance of exciton-exciton inter¬ 
actions, excitation-induced dephasing (EID), and bound 
biexcitons in quantum wells. Recently, the 2DFT spec¬ 
troscopy technique has been applied to semiconductor 
microcavity polaritons m- 

In this paper, we report on 2DFT spectra when both 
lower and upper polariton states are simultaneously ex¬ 
cited. We use 2DFT spectroscopy to differentiate the two 
types of nonlinearities: self and cross-polariton interac¬ 
tions. We perform two-pulse FWM experiments in both 
positive and negative time delay configurations. Polari¬ 
tons inherit, from their photonic component, a light effec¬ 
tive mass that leads to a strong parabolicity in energy- 
momentum dispersion, which is generally neglected in 
bare quantum well excitons. We reveal the role played 
by the energy-momentum dispersion on the nonlinear po¬ 
lariton dynamics, which is usually not involved in 2DFT 
spectroscopy of heavy particles. This paper is organized 
as follows: Section II, describes the sample and the four- 
wave mixing experiment. Section III reports on the ex¬ 
perimental results and a simple third-order perturbative 
analysis, and Section IV is dedicated to a detailed theo¬ 
retical model and numerical simulation using lower- and 
upper-polariton basis Gross-Pitaevskii equations. 

II. EXPERIMENTAL METHOD 

The sample is a high quality III-V GaAs-based micro¬ 
cavity m- A single 8 nm Ino. 04 Gao. 96 As quantum 
well is sandwiched between a pair of GaAs/AlAs dis¬ 
tributed Bragg-reflectors. The Rabi splitting energy at 
zero cavity-exciton detuning (^=0) is Q =3.26 meV. The 
experiments are performed at the cryogenic tempera¬ 
ture of 4 K with several positions on the sample cor¬ 
responding to different exciton-cavity detunings. We use 
a Ti:sapphire laser with a broad spectrum femtosecond 
pulse and 80 MHz repetition rate. The center energy of 
the pulse spectrum is set between the lower and upper 
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polariton energies. We employ four-wave-mixing spec¬ 
troscopy in two-pulses configuration. The k 2 and ki 
pulses arrive on the sample in directions =0.96 
and ^2=0 /im“^ respectively (See Fig. M (a)). The ex¬ 
periments are performed in the low-density regime with 
ki and k 2 pulse intensities of 6.7 x 10^^ photons pulse“^ 
cm“^. The pulses are co-circularly polarized in order to 
avoid the biexciton effect [HI [17] . We detect the FWM 
signal in the direction kpwM = 2^2 — ki = —ki. The 
pulses ki and k 2 arrive on the sample at times and 
tk 2 respectively. A time delay r = — 1^2 between two 

pulses is called positive (negative) when the k 2 {ki) pulse 
arrives before the ki (/C 2 ) pulse. In Fig. [^(b), we show 
schematically the polariton dispersion with the ki and k 2 
pulses and also the k^wu signal. Fig. Be) shows that 
the transmissions of k 2 and ki pulse have different energy 
peaks due to the effect of the polariton energy dispersion. 

The experimental set up is explained in detail in our 
previous work [T8|. With a heterodyne detection tech¬ 
nique [19], we record the electric field of the FWM sig¬ 
nal S'(r, t), where t represents the real evolution time of 
the FWM signal after the incidents of the two pulses. 
Notice that the FWM signal 5'(r, t) is function of two- 
independent time periods, the time delay r and the real 
time evolution t of the signal. We obtain a delay depen¬ 
dent ID FWM amplitude signal spectrum |5'(r, et)| by 
performing a Fourier transformation (FT) with respect 
to t, which is performed by spectral interferometry [6]. 
The spectrometer acts as a FT, converting the real time 
evolution t into the third-order emission energy e^. The 
2D spectrum S{er,€t) is then obtained through the FT 
with respect to both the r and t axes. Here e^- and rep¬ 
resent the absorption and emission energies respectively. 
Before the conversion to the 2DFT spectrum, we apply a 
numerical phase correction [20] using the upper-polariton 
energy as a phase reference (a detailed explanation is 
given in Appendix A). 

III. EXPERIMENTAL RESULTS AND 
PERTURBATIVE ANALYSIS 

In this section, we present experimental results from 
FWM experiments performed with the cavity detuning 
at S = —0.38 meV and analyse them with a conventional 
third-order perturbation theory of nonlinear optics. The 
detail of the third-order perturbation theory is presented 
in Appendix B. 

A. ID FWM spectrum 

In Fig. 1^ (a), we display the amplitude of the ID 
FWM signal |5'(t, et)|, which is the spectrum of the emit¬ 
ted signal as function of the delay time, r, between the 
two pulses. The ID FWM signal spectrum presents two 
main resonances: the lower frequency emission, originat- 


(a) 




FIG. I. (color online) Schematic of FWM configuration and 
pulse sequence (a). Lower- (LP) and upper-polariton (UP) 
energy-momentum dispersion at a slight negative cavity de¬ 
tuning (b). The dashed black lines represent exciton and pho¬ 
ton and energy-momentum dispersion. Transmission spec¬ 
trum of the fe and ki beam (c). eL{k) and eu(k) are respec¬ 
tively lower and upper-polariton enegy at wave-vector k 

ing from lower polariton (LP), and the higher frequency 
emission, from upper polariton (UP). Each one display¬ 
ing a fine structure. Moreover, the FWM emission shows 
a temporal oscillation behaviour with a period of 1.2 ps. 
This period corresponds to the Rabi splitting energy. We 
will show that this oscillation can be understood as a 
quantum beat. All these features are addressed in the 
next sections. 


B. 2DFT spectrum 

The absolute value of the 2DFT spectrum \S{er,et)\ is 
shown in Fig. (b) for negative r and in Fig. (c) for 
positive r. Fig. (b) and (c) are respectively referred 
to as the one-quantum and two-quantum regimes [2]. A 
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Fourier transformation with respect to the time delay |r| 
converts the delay map 5'(|r|,et) into the 2DFT spec¬ 
trum S{er,et). One axis displays the absorption energy 
and the other the emission energy e^. We will analyse 
the 2DFT spectrum in terms of third-order perturbation 
theory with the help of double-sided Feynman diagrams. 


and 

oc -2|r2L|^ 


[i{et — e2L{kFWM) + £L{kFWM)) + 'll] 


[i{er + eL(fci)) + 7 l] 


C. 2DFT: One-quantum regime 

Firstly, we focus on the FWM signal emitted at the 
negative time delay (r < 0). In Fig. (b), we observe 
two diagonal groups, LP-LP and UP-UP, and two off- 
diagonal groups, LP-UP and UP-LP. Inside each peak 
group, fine structures are found. We classify and name 
them the virtual (VB), middle (MB), and normal (NB) 
branches going from lower to higher emission energies. In 
order to analyse the origin of LP-LP and UP-LP groups, 
we introduce double-sided Feynman diagrams (Fig. 
HEIIES]. The double-sided Feynman diagrams repre¬ 
sent Liouville-space pathways of the FWM signals. Each 
diagram corresponds to the third-order perturbative evo¬ 
lution of the system’s density matrix. The FWM signal 
5'(|r|,t) can be calculated, within third-order perturba¬ 
tion, by summing all pathways. In the diagrams, the 
vertical line is the time evolution of the system, with the 
time ordering of the arrival of ki , k 2 pulses and the FWM 
signal emission time. In all pathways, the arrival of the 
first ki pulse creates a coherence between the ground 
state and a single LP ( 0 -LP) or UP (0-UP) state. Both 
the second and third fields come from the pulse k 2 after 
a delay |r| simultaneously, which fixes T = 0 in the di¬ 
agrams. Performing a standard third-order perturbative 
calculation of nonlinear optics [H HO] (See the Appendix 
B), the FWM signal corresponding to the diagram (A) is 
given by. 


has the exactly same form as Here, we treat 

the semiconductor microcavity system as two-oscillators 
(lower and upper-polaritons) weakly coupled to photons 
outside the cavity with the couplings and flu- The 
detailed background of this polariton basis model will 
be discussed in the theoretical model section. Note 
that diagrams (A) and (B) in Fig. include only the 
ground and “single-quantum” state (LP), while diagram 
(C) also includes the “two-quantum” state ( 2 LP). This 
two-quantum state (2LP) is modified by the polariton- 
polariton self-interaction, resulting in the energy of 2LP 
state € 2 L{k) being slightly blue shifted from twice that 
of LP state 2eL{k) {i.e., € 2 L{k) 7 ^ 2eL{k)). It is worth 
noting, that if the lower-polariton self-interaction were 
absent, the relation e 2 L{k) = 2eL{k) would hold and the 
sum would be zero. This is an intu¬ 

itive consequence of the fact that no FWM signal appears 
in a linear system [23]. The same description and path¬ 
ways are applied to the UP-UP resonance, considering, 
in this case, only the upper polariton coherence and the 
upper-polariton self-interaction. 

Similar to the LP-LP group, we present the Liouville- 
space pathway processes of UP-LP FWM signals in terms 
of double-sided Feynman diagrams in Fig. @(D)-(G). 
The signals associated with these diagrams are given by. 


S'(^’^)(e^,et) oc 


__ 

[i{et — eLikpwM)) + "fL][i{^T + + ^u] 


where eL{k){eu{k)) is the energy of the lower (upper)- 
polariton, ^^(Ut/) represents the coupling constant be¬ 
tween the lower (upper)-polariton and the photon out¬ 
side the cavity, and is the dephasing rate of lower 

(upper)-polariton. In this pathway, during |r| and t, the 
system evolves keeping a coherence between its ground 
state and a single lower-polariton state ( 0 -LP and LP-0). 

In the “one-quantum regime”, the first order evolution, 
during time |r|, is always a coherence between the ground 
state and single lower or upper-polariton state (0-LP or 
0-UP). The Fourier transformation of 5'^^^(|r|,t) reads, 

oc ^.- 

[i{et - eLikpwM)) + + 7^] 

Similarly, the contributions from diagrams (B) and (C) 
lead to 


and 

OC —' 


[i{et - eLui^EWM) + ^u{^fwm)) + 11] 


S^^\er,et) oc 




[i{et — eLikpwM)) + 7L][^(er + ^L{ki)) + 7l] 


[i{er + + lu] 

Along pathways (F) and (G) the FWM emission origi¬ 
nates from the coherence between the two-quantum state 
(LP+UP) and the single-quantum state (UP). Similar 
to the pathways (A)-(C), the energy of UP-LP state 
is shifted due to the lower and upper-polariton cross¬ 
interaction, €Lu{k) 7 ^ €L{k) €u{k). Again, if the 

' lower and upper-polariton cross-interaction does not ex¬ 
ist, eLu{k) = eL(^) + ^u{k) holds and the summation of 
the pathways (D)-(G) cancels. This leads to the disap¬ 
pearance of the off-diagonal peaks. We can draw simi¬ 
lar diagrams for LP-UP groups (not shown) and calcu¬ 
late perturbatively the FWM signals. The plot of cal- 
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FIG. 2. (color online) Experimental amplitude of FWM spec¬ 
trum as a function of emission energy and k 2 -ki pulse delay 
r: et)| (a). Amplitude of 2DFT spectrum |5'(er,et)| for 

1-quantum (b) and 2-quantum region (c). Diagonal dashed 
lines represent Cr = ct for (b) and Cr = 2et for (c). Horizontal 
and vertical dashed lines respectively represent different ab¬ 
sorption and emission energies. ^ virtual branch 

(VB), which is explained in Section IV. c. Colour scales are 
normalized by the maximum and minimum of the amplitude. 
The same normalization and colour bar are used in all figures 
of the article. 



(b) 2-Quantum 

E'ii-k,) Em Eii-k,) E'u{-k,)Eu{0)Eu{-k,) 



2Eu{0) 


£l(0)+£u{0) 


2El{0) 


FIG. 3. (color online) Plot of the amplitude of 2DFT spec¬ 
trum |AS'(er, et)| based on third-order perturbative calculations 
for the (a) 1-quantum and (b) 2-quantum regimes. 


double-sided Feynman diagram analysis elucidates the 
origin of the amplitude oscillation in the FWM emission 
along the delay of the two pulses (Fig. (a)). This 
amplitude oscillation can be understood as a quantum 
beat: an interference of the pathways (A)-(C) and (D)- 
(G). During the delay r, in the pathways (A)-(C) the 
phase evolves as while in the pathways (D)-(G) 

it evolves as Thus, the amplitude of the FWM 

signal of the delay map (Fig. (a)) shows a beat fre¬ 
quency which corresponds to the Rabi splitting energy: 
eu — eL — 3.26 meV. 


culated 2DFT spectrum including all pathways in the 
one-quantum regime is shown Fig. [^(a). 

As we discussed above, the diagonal peaks arise 
from the polariton self-interactions while the off-diagonal 
peaks arise from the cross-interactions. Moreover, the 


It is important to note that the lower and upper- 
polariton branches have energy-momentum dispersions: 
eL(fc) ~ CLfi + and ejj{k) ~ eu,o + (S®® 

Fig. 0 (b)), where tul and mu are the mass of the lower 
and upper-polaritons respectively. As €L(u){kFWM) = 
€L(u){~ki) = ^L(u){ki), the LP-LP (UP-UP) peaks are 









































5 


LP-LP 

(A) 


(B) 


sidered [24]. For doing this, we employ non-perturbative 
numerical simulations in Section IV. 



0 




FIG. 5. Double sided Feynman diagrams that represent per¬ 
turbative Liouville-space pathways of 2-quantum 2DFT spec¬ 
trum. Diagrams (H-I) and (J-M) respectively represent the 
normal branches of 2LP-LP and LPUP-LP group (Fig. [^(b)). 
Polariton-polariton interactions are present in all pathways. 
ql and gx respectively represent self and cross-interaction of 
polaritons. 


FIG. 4. Double sided Feynman diagrams that represent per¬ 
turbative Liouville-space pathways of 1-quantum 2D FT spec¬ 
trum. Diagrams (A-G) and (D-CD respectively represent the 
LP-LP and UP-LP peaks (Fig. [^ (b)). Polariton-polariton 
interactions are introduced in pathways (G), (F), and (G), 
where double k 2 pulses excite two-polariton state, ql and gx 
respectively represent self and cross-interaction of polaritons. 

absorbed and emitted at the same energy in the third- 
order perturbation theory. These peaks presented in 
Fig. [^a) correspond to the normal branches (NB) of 
the experimental 2D spectra in Fig. (b). We notice 
that the third-order perturbative model reproduces only 
the normal branches, which are resonant to the polari- 
ton energy-momentum dispersion, and no fine structure 
appears inside each group. This is because polariton- 
polariton interaction is considered only as a level-shift 
of the eigen state energy in the third-order perturbation 
theory. In order to give rise to fine energy structures, 
the superposition between different momentum states in¬ 
duced by polariton-polariton interaction needs to be con- 


D. 2DFT: Two-quantum regime 

We concentrate now on the FWM signal emitted at 
positive time delays (r > 0), where the k 2 pulse arrives 
first. Since the k 2 pulse acts as two degenerate pulses, 
this pulse creates a coherence between the ground state 
and the two-quantum state. Thus we call this time de¬ 
lay configuration “two-quantum regime”. Double-sided 
Feynman diagrams corresponding to 2LP-LP and LPUP- 
LP groups are presented in Fig. In all pathways, the 
system evolves in two-quantum coherence (0-2LP or 0- 
LPUP) during time r, before the arrival of the second 
pulse ki. For example, the FWM signal corresponding 
to the diagram (H) can be written as 
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For which the Fourier transformation gives, 
S^^\er,et) oc 


[i{et - eLikpwM)) + 7l] 

1 


[i(er - e2L(^2)) + 27 l] 


Similarly, 
S^^\er,€t) 
oc -2|r2L|^ 


[i{€t — €2L{kFWM) + ^LikpWM)) + 7i] 

1 


[i{er - €2L{k2)) + 27 l] 

For LPUP-LP groups, the FWM contributions read, 
S^^\er,et)^\nL\^\nu\^- ^ 


[i{et — ^hikpwM)) + 7l] 

1 


[i{^r - ^L{k 2 ) - eu{k 2 )) + 7 l + 7c/] 


Similarly, 
S^^\er,et) 
oc — 


[i(et - e2L{kFWM) + ^l{^fwm)) + 7 l] 

1 


[/(Cr - eL(/i: 2 ) - ^u{k 2 )) + 7l + 7c/]' 


S^^\er,et) oc jf^Lpl^c/l^ 


[i(et - eLikpwM)) + 7l] 

1 


[/(er - eL{k 2 ) - eu{k 2 )) + 7l + 7c/] ’ 


and 

(e^, et) 

^ “I^Lpl^C/P 


[/(er - ^L{k 2 ) - eu{k 2 )) + 7l + 7c/] 

Here, we do not repeat the same discussion for the LPUP- 
UP and the UP-UP groups. We plot the calculated 2DFT 
spectrum of the two-quantum contribution from all path¬ 
ways in Fig. [^(b). Energy shifts, originating from self 
and cross-interactions e 2 L{k) ^ ‘^CLik) and eLu{k) ^ 
^hik) -\-eu{k)^ are necessary for the appearance of the on 
and off-diagonal peaks respectively. Again, the diagonal 
and off-diagonal groups are associated with the self and 
cross-interactions respectively. Notice that the diagonal 
line in the figure (e^- = 2et) is defined with an absorption 
energy that is twice the emission energy, this is charac¬ 
teristic of a “two-quantum regime”. In Fig. (b), the 
dashed diagonal line does not pass through the normal 
branches (NB) of 2LP-LP and 2UP-UP peaks. This is a 
consequence of the polaritons energy-momentum disper¬ 
sion: eL{U){^FWM) = ^L{U){ki) > eL{U){^ 2 ) = ^L{U){^) 
(See Fig. 3(b)). 


IV. THEORETICAL MODEL 

Our starting point is the bosonic exciton Hamiltonian 
including the exciton-exciton interaction and the exciton- 
photon coupling, 

Hex ~ Huyi T Hxint' (1) 

The linear term Hun and nonlinear exciton-exciton in¬ 
teraction term Hxint are respectively given by 


7' 


Hiin — / dr 


fc 2^2 ^ ^ ^ 

- -7 ,— )4’x + i’lF - -7 — )i’c 


2mr 


2mc 


+ ^(^17 + '^I'^i’x) + ^(7V’6 + 


and 


Hxir 




( 2 ) 


(3) 


[i{et — e2L{kFWM) + ^hi^FWivt)) + 1l] 


ipx and ipc are the exciton and cavity pho¬ 

ton annihilation (creation) operators respectively. The 
Rabi splitting between exciton and cavity photon is rep¬ 
resented by O. flqm is a quasi-mode Rabi splitting which 
is the coupling of photons between the inside and outside 
of the cavity [25]. -05 represents a classical photon field 
outside the cavity. Cx and Cc are respectively the exciton 
and photon eigenenergies. The exciton-exciton interac¬ 
tion constant is given by go [26]. Now, we introduce the 
polariton bases tpL and defined as, 

^x\^(X -C\(^L 

to rewrite the Hamiltonian Eq. [^ X and C are re¬ 
spectively excitonic (photonic) and photonic (excitonic) 
Hopfield coefficients of lower (upper) polaritons. They 
are chosen to be. 


and 



C = -x F 1- 




which diagonalizes a non-interacting exciton-photon 
Hamiltonian at k=0. We rewrite Hun in terms of po¬ 
lariton basis under a parabolic approximation of polari¬ 
ton energy-momentum dispersion: eL(u)(/c) 6^,0 + 

Hamiltonian of the linear part in polari¬ 
ton basis reads, 


Hlin ^ Hi 


lin 


= dr 


/ 




+ ^^^6 + 


lU 


( 4 ) 
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and et/^o are respectively the energies of lower and 
upper-polariton at zero momentum written as, 

eL,o = ^ (2e, + (5 - 

and 

The polariton quasi-mode Rabi splittings are defined as 
VLl = CVtqm and Vtu = X^tq^. The polariton mass is 
given by l/mnu) = \X\'^/m^(c) + \C?/mc(x)- Now, the 
linear Hamiltonian Eq. is formally the same as that 
of two oscillators, which are weakly coupled to photonic 
fields. Actually, this is why the double-sided Feynman 
diagram analysis, which is normally used in a weak- 
coupling system between oscillators and photons, can 
be applied to this system in which exciton and photon 
are strongly coupled inside a microcavity. We introduce 
polariton-polariton interactions as, 

Hpint = + ]j^9U'ti’\j'^\j'^U'^U- 

■ (5) 

This is a simple effective Hamiltonian that can give 
self- and cross-interactions between lower and upper- 
polaritons. Using the Hopfield coefficients, we set ql = 
go\X\‘^ , gu = go\C\^, and gx = 2go\X\'^\C\'^. The rela- 
tion between this Hamiltonian and exciton-exciton inter¬ 
action Hamiltonian {Hxint) is discussed in Appendix C. 
Only when the kinetic term (energy-momentum disper¬ 
sion) in Eq. is neglected, can we derive the effective 
polariton-polariton Hamiltonian {Hpint) from exciton- 
exciton interaction Hamiltonian {Hxint) using a pertur¬ 
bation theory in the low density regime. 

In order to understand the experimental observations 
in detail, we perform non-perturbative numerical simu¬ 
lations. A polariton basis total Hamiltonian is, 

Hpoi = Hi^^ H- Hpint- (b) 


With the aid of the Heisenberg equations of motion 
Hpol] and the mean-field approxi¬ 
mation the equations of motion of lower and upper 
polariton wavefunctions are simplified to non-equilibrium 
lower-upper polariton Gross-Pitaevskii equations: 


ih'ipL = (el.o - +flL|V’Lp 

2mL 

+9x\^U? - + ^/ext 

h? 

ihipu = (et7,o - X- 

Zmu 


(7) 


( 8 ) 


where ipL{u) = {i^L{u)) is lower (upper) polariton wave 
function. The polariton decay rate is given by 7 l(c/) = 
\^?‘lx{c) + |C'p7c(cc), where ^x and 7c are chosen to be 
the same (0.33 meV). The k 2 and ki pulse excitations 
are represented by an external photon field /ext(= '^b)- 
The constant gL{u) represents a self-interaction of lower 
(upper) polaritons, while gx is a cross-interaction con¬ 
stant between the lower and upper-polariton. In analogy 
with nonlinear optics, gL{u) und gx can be called the 
self-phase modulation (SPM) and the cross-phase mod¬ 
ulation (XPM) terms respectively. Similar to nonlinear 
optics, the XPM term gx is twice as strong as the SPM 
term gL{u) EH]- lu the simulation, the exciton-exciton 
interaction constant go is set as 2 meV/no, where no is a 
normalized density. The excitation |0*^/ext is a Gaus¬ 
sian pulse with a peak intensity of 0.5y^ and a pulse du¬ 
ration of 250 fs. This model describes coherent processes 
and includes mean field interaction. Using this simplified 
model, we can directly investigate the contribution of self 
and cross-polariton interactions through the three nonlin¬ 
ear interaction constants: gL^ gu^ and gx- This is a key 
advantage of this model over the exciton-photon basis 
(local mode basis) Gross-Pitaevskii equations (See Ap¬ 
pendix G). In the local mode basis, the only nonlinear in¬ 
teraction constant represents the exciton-exciton interac¬ 
tion, thus we cannot deal with self and cross-interactions 
of lower and upper polaritons independently. Simulated 
FWM spectra are presented in Fig. All numerical sim¬ 
ulations are performed in one-dimensional space. Simi¬ 
lar to experimental observations (Fig. [^(b) and (c)), the 
simulated 2DFT spectra (Fig. [^(b) and (c)) clearly show 
fine structures inside the four peak groups. 


A. Different interaction contributions 

In order to obtain better insight into the importance 
of the interactions, two different sets of spectra are cal¬ 
culated and plotted in Fig. Firstly, Fig. (a) shows 
the 2D spectra when considering only self-interactions 
9l{u)' As expected, we find two main LP-LP (2LP- 
LP) and UP-UP (2UP-UP) groups along the diagonal 
line in one (two)-quantum 2DFT spectrum. The LP- 
LP (2LP-LP) and UP-UP (2UP-UP) groups originate 
from lower-lower and upper-upper self-interactions re¬ 
spectively. On the other hand, the 2D spectra, including 
only the lower-upper cross-interaction gx (Fig. [!|(b)), 
shows only the off-diagonal groups UP-LP (LPUP-LP) 
and LP-UP (LPUP-UP) in one (two)-quantum 2D spec¬ 
trum. Only when we include both the self gL{u) cross 
gx interactions, is the observed experimental 2D spectra 
are reproduced (Fig. [^. Notice that, since the cross¬ 
interaction constant gx is twice as strong as the self¬ 
interaction constant gL{u)i around zero cavity detuning, 
the off-diagonal peaks are brighter than diagonal ones, 
corroborating the experimental results. 






(a) Delaymap 

£'lW £l{0) £L{-k,) 


(a) SPM 


(b) XPM 


£u{-k,) £u{0)£u{-k,) 



- 3 - 2-10123 

(b) 1-Quantum 

£'L(-k,) £l{ 0) £L{-k,) £'u(-k,)£u{0)£u{-k,) 



n- 1 -r 

- 2-10 1 
emission energy a (meV) 


£u{k,) 


£L{k,) 


2£u{0) 


2£l{0) 


FIG. 6. Simulated amplitude of FWM spectrum as a fuuctiou 
of emissiou euergy aud fe-Zci pulse delay r: et)| (a). 

Simulated amplitude of 2D FWM spectrum |5'(er, et)| for the 
(b) 1-quautum aud (c) 2-quautum regious. 


Delaymap 


Delaymap 



emission energy a (meV) 


emission energy a (meV) 


FIG. 7. Simulated delaymap |5'(t, et)| aud 2D FT spectrum 
|AS'(er,et)| for 1-quautum aud 2-quautum regime. Simulatious 
are performed iucludiug ouly self-iuteractiou (or self phase 
modulatiou: SPM) gx = 0 (a) aud for ouly cross-iuteractiou 
(or cross phase modulatiou: XPM) gL(u) = 0 (b). 

laritons energy-momentum dispersion, associated with a 
light polariton mass. The idea is schematically shown 
in Fig. 1^ (one-quantum regime) and (two-quantum 
regime) as energy diagrams. In this section, we focus 
on the fine structure in the one-quantum regime how¬ 
ever, the discussion is completely the same for the two- 
quantum regime. 


B. Fine structures inside peak groups 


In the 2DFT spectra, we can find a fine structure (the 
normal (NB), middle (MB), and virtual branch (VB)) 
inside each peak group. They can be classified by the 
emission energies. For example, in the LP-LP group, the 
emission energies of NB, MB and VB respectively cor- 
respond to + ^^kl, eLfi and eLfi - ^^k^. As 
mentioned above, the fine structure is related to the po- 
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eL(k) 



£L(k)+£L(k) 

SL(k) 


(II) VB 



£L(k)+£L(k) 

£L(k) 



(III) MB 

£L(k) 



UP-LP 



time-ordering 



First Second FWM 

pulse pulse signal 




(III) MB 


UP-LP 







FIG. 8. Schematic energy diagrams representing the origin 
of the fine structures in one-quantum 2D spectrum. Three 
solid arrows represent two degenerate k 2 pulses and ki pulse. 
The dashed arrows are FWM emissions. For NB schematics, 
the indices (A)-(G) correspond to those of the double-sided 
Feynman diagrams in Fig. 

Firstly, we discuss the process associated with the nor¬ 
mal branches. As is shown in Fig. gl) and (IV), NB 
emission is an on-branch FWM emission. Therefore, 
its emission energy is eL(u){kFWM) = ^L{U){-ki) = 
^L(( 7 ),Q+ 2 m^(^) both diagonal and off-diagonal peak 
groups. In the third-order perturbative calculation of 
FWM signal in Section III, only this branch appears (See 
Fig. §, because the FWM signal evolve with the eigenen- 
ergy of polaritons, which corresponds to (See 

Appendix B). 





FIG. 9. Same schematic as Fig. [^for 2-quantum 2D spec¬ 
trum. For normal branch (NB) schematics, the indices (H)- 
(M) correspond to those of the double-sided Feynman dia¬ 
grams in Fig. 


Secondly, let us consider the virtual branches. The ap¬ 
pearance of the VB can be understood in a framework 
of off-resonant scattering induced by polariton-polariton 
interactions m- The fine structure results from mo- 
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mentum and energy conservation in the FWM process. 
As the energy diagrams of the VB, in Fig. (II) and 
(V) show, if the ki and k 2 pulses are on-resonant exci¬ 
tation, energy and momentum conservation demand VB 
emission. Using energy conservation, for LP-LP group 
(See Fig. 8 (II)), the energy of VB emission is given by 
= 2eL(A:2) - = 6^,0 - 2^^1- Si^ce 

NB emission energy is 6^,0 + separation be- 

tween NB and VB emission is which is twice 

the energy difference between the ^2 and ki transmis¬ 
sion peak of the lower polariton. In the experiment, 
the energy separation between NB and VB emission is 
0.58 meV (See Fig. (b) and (c)). This energy differ¬ 
ence corresponds to the twice of 0.29 meV, which is the 
energy difference of the ki and k 2 transmission peak of 
the lower polariton presented in Fig. (c). Addition¬ 
ally, this relation holds for both for negative and pos¬ 
itive cavity detuning (see Sect. IV C). Therefore, we 
can safely say that VB peak in LP-LP group originates 
from this process. The VB energy-momentum disper¬ 
sion e^(/c) is a mirror image of NB energy-momentum 
dispersion eL(k). We note that in the weak intensity 
excitation regime, the dispersion is still parabolic and 
we do not consider linearisation of the dispersion due 
to the Bogoliubov transformation [ 6 ]. In UP-LP group 
(See Fig. |^(V)), similar to the LP-LP group, from mo¬ 
mentum and energy conservation VB emission energy is 

^Lui^FWAd) = ^L{k2) + eu{k2) — eu{ki) = eL,o - 
Notice that, in Fig. |^(V), the higher energy-momentum 
dispersion is the sum of two different energy-momentum 
dispersions eL{k) and €u{k). While the VB of the UP-LP 
group is visible in the simulation. Fig. experimentally 
we cannot find the VB in the UP-LP group (See Fig. 
[^. The reason for this is not clear. We can apply the 
same discussion to the LP-UP and the UP-UP groups, 
but experimentally the VB cannot be found in the 2DFT 
spectra of these groups. 

Now, let us focus on the middle branches. In the UP- 
LP group of the experimental 2DFT spectra, a strong 
peak is found next to NB. The emission energy of this 
peak corresponds to eL{k 2 ) = 6 ^, 0 , which is located in 
the middle of NB and VB. Thus we name this peak MB. 
In Fig. we can identify a weak MB in LP-LP. Addition¬ 
ally, the experimental peak in UP-UP is also considered 
as a MB because the emission energy of the UP-UP peak 
is eu{k 2 ) = eu,o- On the other hand, the numerically sim¬ 
ulated 2DFT spectra does not include MBs, which makes 
the interpretation of MB difficult. Finally, schematically, 
we associated the MB to the processes described in Fig. 
ini) and (VI). These are the processes where energy 
and momentum conservation are satisfied and one of the 
degenerate fields of the second pulse is not resonant to 
the energy-momentum dispersion. To understand the de¬ 
tailed mechanism of the MB, a further investigation and 
model are necessary. For instance, taking into account 


effects such as, excitation induced dephasing (FID) and 
relaxation of the upper-polariton into exciton reservoirs. 


C. FWM spectra at S=-2 meV and S=2.7 meV 

To obtain more insight into the origin of the fine struc¬ 
ture energy, we discuss the results obtained at both nega¬ 
tive S=-2.0 meV and positive 6=2.7 meV cavity detuning. 
In Fig. we display the amplitude of the 2DFT sig¬ 
nal \S{er,€t)\. We observe the same features presented 
in Fig. [^(b) and (c), however the emission energy sepa¬ 
ration between the NB and the VB varies depending on 
the cavity detuning due to the change in the polaritons 
energy-momentum dispersion. 


Cavity detuning -2 meV Cavity detuning 2.7 meV 

Delaymap Delaymap 




1 Quantum 



1.486 1.488 1.490 


2 Quantum 



1.486 1.488 1.490 

emission energy a (eV) 


FIG. 10. Experimental FWM delaymap |*S'(t, et)| and 1- and 
2-quantum 2D FT spectrum |AS'(er,et)| for negative (6 = —2 
meV) and for positive {6 = 2.7 meV) cavity detuning. 


By detuning the cavity photon energy far below the 
exciton energy S=-2.0 meV, the LP becomes photon-like 
and acquires a lighter mass than the UP. Due to the 
lighter mass of the LP at S=-2.0 meV, the NB and VB 
in the LP-LP group are more separated than in the ex¬ 
perimental results at the cavity detuning 6 = —0.38 meV 
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(Fig. I (a)). 

On the other hand, positive detuning 5 = 2.1 meV re¬ 
sults in a flatter energy-momentum dispersion for the LP. 
The resulting heavier mass makes the fine structures in 
the LP-LP group almost degenerate. We also find three 
fine structure peaks in UP-UP (2UP-UP) group. This 
is due to the lighter mass of the UP with positive cavity 
detuning. However, the brightest MB peak cannot be de¬ 
scribed within our simple model. It might be related to 
dynamic redshift of the UP mean field energy [6l [18] . We 
also observe weaker spectral features in the one- and two- 
quantum 2D spectra, which could originate from higher 
order non-linear processes that were not completely elim¬ 
inated in our heterodyne detection scheme. This will be 
the subject of further investigations to confirm their ori¬ 
gin. 

The simulations for negative {S = —2 meV) and pos¬ 
itive {6 = 2.7 meV) cavity detuning are shown in Fig. 
TTj As expected, for the negative detuning the energy 


separation of fine structures (NB and VB) in the LP- 
LP group increases because the LP dispersion becomes 
photon-like. Conversely, for positive detuning, we cannot 
distinguish the fine structure inside the LP-LP group due 
to the flat dispersion of the LP (exciton-like). The fine 
structures in UP-UP group shows the inverse detuning 
dependence as the LP-LP group: a small separation for 
S = —2 meV and a large separation for S = 2.7 meV. 
Notice that the energy separations between NB and VB 
in UP-LP (LPUP-LP) and LP-UP (LPUP-UP) groups 
are the average of those of LP-LP (2LP-LP) and UP-UP 
(2LP-LP) groups. 


D. Comparison with experiments 


Cavity detuning -2 meV Cavity detuning 2.7 meV 

Delaymap Delaymap 




1 Quantum 



2 Quantum 



-2 -1 0 1 2 3 4 5 


emission energy a (meV) 


FIG. 11. Simulated FWM delaymap |*S'(t, et)| and 1 and 2- 
quantum 2D FT spectrum |5'(er,et)| for negative {S = —2 
meV) and for positive {S = 2.7 meV) cavity detuning. 


Our numerical model qualitatively explains the ap¬ 
pearance of LP-LP (2LP-LP) and UPLP (LPUP-LP) 
peak groups and the fine structures inside them. In prin¬ 
ciple, as is shown in the simulations, the fine structure 
should appear in the group LP-UP (2LP-UP) in the ex¬ 
perimental spectrum. However, experimentally there is 
a strong amplitude asymmetry between UP-LP (LPUP- 
LP) and LP-UP (LPUP-UP) group. Compared with UP- 
LP (LPUP-LP) group, the LP-UP (LPUP-UP) group 
is very weak and we cannot resolve detailed structures. 
This type of asymmetry in the off-diagonal peaks is ob¬ 
served also in Ref. m and in 2DFT spectra of bare 
quantum well [T4l[22| experiments. Additionally, the fine 
structure in UP-UP group is more complicated than what 
is computed within our simple model. The theoretically 
predicted off-diagonal VBs do not appear in the experi¬ 
ment. On the other hand, MBs, which are not visible in 
the simulation, appear in the experiment. 

Although a microscopic model completely explaining 
the experiments is still lacking, in term of asymmetric 
amplitudes of off-diagonal peaks, we can attempt to fit 


the observed spectra by dealing with the polariton Gross- 
Pitaevskii equations in a more phenomenological way. 
Here, for the fitting between the experiment of the cavity 
detuning S = —0.38 meV and simulation, we consider the 
strength of the interaction constants of polaritons as free 
fitting parameters: 

ihi^L = (cl.O - 

2mL 

+fl[/Ll^i7p-*Y)V’L + ^/ext (9) 

fl2 

ihipu = {eu,o - - - 9u\'^uf 

2mu 

+9'lu\'<Pl?+ ( 10 ) 

The above equations are formally same as Eq. 
However, the interaction constants (^j^, g'jj^ and 

g'jj]^) are fitting parameters and no more connected to 

the exciton-exciton interaction Hamiltonian Hint- This 
means that g'j^jj is not necessary equal to g'jjj^. This is 
the advantage of using lower- and upper-polariton basis 
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equations compared with a conventional exciton-photon 
Gross-Pitaevskii equations, where the amplitude of LP- 
UP group is always same as that of UP-LP groups. We 
set the interaction constants as g'j^ : g[j : g'j^jj : g'jjj^ = 1 : 
0.6 : 0.3 : 1.6 (meV/no), which correspond to the ratio of 
the integrated amplitudes of peaks groups LP-LP, UP- 
UP, LP-UP and UP-LP in Fig. [^b). The other param¬ 
eters are the same as those in the previous section. The 
simulated 2DFT spectra with these parameters are shown 
in Fig. We find that our phenomenological model re¬ 
produces more closely the asymmetric intensities of off- 
diagonal peaks. The relative strength of the interaction 
constants is believed to reflect neglected contributions 
such as a frequency dependent non-Markovian nature of 
exciton-exciton interaction (exciton-exciton correlation) 
[151 [ini ISO] and photon-assisted exchange interaction be¬ 
tween polaritons that reinforces the repulsive interaction 
among the lower polaritons but that weakens the repul¬ 
sive interaction among the upper polaritons Eiisi]. 

V. CONCLUSION 

In summary, 2D Fourier transformation spectroscopy 
is performed to investigate polariton-polariton interac¬ 
tions in semiconductor microcavities. The experimen¬ 
tal 2D optical spectra demonstrate the existence of 
lower-upper polaritons cross-interaction and of lower (up¬ 
per) polariton self-interaction, which originate from the 
Coulomb and exchange interactions between the fermion 
constituents of the exciton-polariton. Furthermore, an 
asymmetry of the coupling between the upper and lower 
polaritons is clearly evidenced in these spectra and indi¬ 
cates complex many-body effects such as exciton-exciton 
correlation and photon-assisted exchange scattering be¬ 
tween carriers constituting the exciton-polaritons. In ad¬ 
dition, a fine structure in the emission energy is iden¬ 
tified as resulting from the polariton energy-momentum 
dispersion and the optical non-linearity of the third or¬ 
der. This work opens the way for a quantitative study 
of many-body effects on composite bosons based on two- 
dimensional Fourier transform spectroscopy. 
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APPENDIX A: NUMERICAL PHASE 
CORRECTION 

Due to the lack of phase stabilization in our experi¬ 
mental setup, we apply a numerical phase correction 


(a) e'Lik,) Sl{o) Stik,) e'u{o)au{o) Su{k,) 



(b) 1-Quantum 

e'Lik,) Siio) SLik,) e'u{o)eu{o) Suik,) 
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(c) 2-Quantum 
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FIG. 12. Simulation based on interaction constants as free 
parameters. The cavity detuning is (5 = —0.38 meV. The 
interaction constants set as g^ ■ g'u ' g'hu • 9ul — 1 • 0.6 : 0.3 : 
1.6 (meV/no). Amplitude of FWM spectrum as a function 
of emission energy and k2-ki pulse delay r: |aS'(t, et)| (a). 
Simulated amplitude of 2D FWM spectrum |*S'(er, et)| for the 
(b) 1-quantum and (c) 2-quantum regions. 

process to the 2D spectrum m . Phase corrected 
2D spectrum Scor can be obtained through the 
Fourier transformation of S'cor(ct,T), where 5'cor(ct,T) = 
S'(et,r)exp(|{ecorT - arg[S'(ecor, t)])). Here ecor repre- 
sents a phase correction frequency. In the article, we 
chose the upper polariton energy as the phase correction 
energy. In the case of Fig. where cavity detuning is 


















13 


S = —0.38, ecor is equal to 1.4883eV. 

In Fig. we present experimental and simulated 
two-quantum 2D spectra with different phase correction 
energies. Fig. 13 (b-c) and (d-e) are respectively two- 
quantum 2DFT spectra using the NB in upper-polariton 
and VB in lower-polariton as phase correction energies. 
Both simulated and experimental spectra indicate that 
phase correction process mainly shifts the axis and 
affects the amplitude of 2D spectrum, but it does not 
change the fine structures of the peaks. 


(a) without phase correction 




emission energy a (meV) 



FIG. 13. Simulated two-quantum 2DFT spectra |5'(er,et)| 
without phase correction (a), which is same spectrum used 
in Fig. 1^ (c). Simulated (b,d) and experimental (c,e) 
two-quantum 2DFT spectra |AS'(er,et)| for choosing upper- 
polariton (b-c) and lower-polariton (d-e) as phase correction 
energies. Arrows represent the phase correction energies Ccor- 


Appendix B: Third-order perturbation theory 

In this appendix, we briefiy explain the calculation of the 
FWM signal based on third-order perturbation theory. 
Firstly, we introduce an ideal three-levels exciton system 
interacting with classical electric fields. \g) represents a 
ground state. A first |e) and second excited states \m) are 
respectively one-exciton and two-excitons states. Now, a 
density matrix can be constructed as following: 

( Pgg Peg Pmg \ 

Peg Pee Pme 1 

Pmg Pme Pmm j 

The system’s Hamiltonian R is composed of an eigen 
Hamiltonian i^o and an exciton-electric field coupling g: 
H = Ho ^ g. Ho and g are given by, 

/ e, 0 0 \ 

Ho=\ 0 Ce 0 ( 12 ) 

\ 0 0 Cm J 

and 

/ 0 df;* 0 \ 

li = 0 y/2VtE^ . (13) 

V 0 ^/ 2 ^ 1 *F; 0 / 

e^, Ce, and are respectively the ground, first, and sec¬ 
ond excited state energies. The time evolution of the 
density matrix p is determined by the Liouville-von Neu¬ 
mann equation ihp = [H^p] = Hp — pH. Firstly, without 
considering the momentum degree of freedom, we can 
obtain 6 coupled equations of motion for the density ma¬ 
trix elements. Then, we perturbatively expand the den¬ 
sity matrix in terms of the orders of the incident electric 
fields as 


_ (0) , (1) , (2) , (3) , 

Pij — Pij ^ Pij ^ Pij + Pij +••• 


(14) 


The first, second and third-order density matrices are 
given by 

= (ee - eg - h)p^ej + 


ihpfj = -in*Ep(^> - nE*p(^J) 

ihp^^} = -iVep^e + {n*EpW* - nE*pW) 

ifipml = (em - eg - n)/3mg + V2n*Epi]^ (15) 

= (ee - eg - i7)pig^ 

+n*E{pg^ - pfj) + V2^E*pgl 

ihp^l = {em - ee - il)ptl + V2n*Ep^^J - ^E*pg\. 

Fe and 7 respectively represent phenomenological decay 
rates of the population of the first excited state and 
polarization. Now we introduce the degree of freedom 
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of momentum and extend the density matrix into 9x9 
taking into account three different momentum channels: 
^ 2 (= 0), ^ 1 , and kFWM{= —ki). We refer to the three 
channels respectively as “pump”, “trigger”, and “idler”. 
Since a superposition state between two different momen¬ 
tum channels does not appear within the following third- 
order perturbation calculation, the 9x9 density matrix is 
block diagonalized into three 3x3 matrices (pump, trig¬ 
ger, and idler). We substitute the density matrix. 


(n) _ (n) 

Pij ~ Pij,p 


(n) iki-x 
■ Pij,t^ 




(16) 


and electric fields. 


E = Ep{t) + Et{t)e^^^^. (17) 

into Eq. and select combinations that satisfy mo¬ 
mentum conservation (phase-matching condition). The 
two pulses Ep and Ef respectively correspond to k 2 and 
ki pulses in the main text. These combinations dif¬ 
fer between the one-quantum (negative delay) and two- 
quantum regime (positive delay). Firstly, in the one- 
quantum regime, there are three pathways A, B, and C: 


path A 

= (4 - 4 “ n)4b + 

= -^*Epp^^> (18) 

^^pfgf = (4 - 4 “ 4)4g!4 + ^*EpP^gg,i 

path B 

= (4 - 4 - 4)4b + ^*^tpflp 

ihpfl^ = -iT,p% + (19) 

= (4 - 4 - 4)4?,f - 

path C 

= (4 - 4 “ 4)p4!t +^*Etpflp 

ihpfi = -iVepfi +( 20 ) 

= (4 - 4 - n)p4:f + V2n*Epp^^l 


Here, 4(m)’ 4(m) respectively represent eigen 

energies of the pump, trigger and idler momentum: 
4(„) = ee(„)(fc 2 ), 4(m) = ee(m)(fci), and = 

^e{m){kFWM)- Secondly, in the two-quantum regime 
(positive delay), two pathways H and I exists: 


path H 

= (4 - 4 “ 4)p44 + ^*Eppfip 
ifip^^ip = (4 - 4 - 4)4Lp + V2n*Eppg\p (21) 
= (4 - 4 “ ^^)pfgf + 4211^(44,P 


path I 

= (4 - 4 “ 4)44 + ^*^ppflp 
if^Pml,p = (4 - 4 “ 4)4s,p + ^^*EpP^elp (22) 
= (4 - 4 - 4)4e,i - p 

The labels A-I correspond to those in the main text. 
These pathways are represented by the double-sided 
Feynman diagrams presented in Fig. and Finally, 
with the third-order polarization density matrices, the 
third-order polarization is calculated as 

+ c.c. = Tr[pf V] 

= npf\ + V2np% + c.c., (23) 

Now, we directly integrate Eq. [TSp^ and calculate the 
third-order polarization density matrices. In general, the 
two pulses are written as 


Ep{t) = 

(24) 

Et{t) = 

(25) 


Ep{t) and Et{t) are respectively pulse envelopes. €pu(tr) 
and tp(^t) are the center energy and arrival time of the 
pulse. For an analytical integration, we assume the two 
pulses have delta function envelopes: 

Ep{t) = EpS{t — tp) and Et{t) = E^S{t — E). (26) 


Let us consider one-quantum regime and the path A. We 
can directly integrate Eq. [^as following, 

J — oo 

— 4o* fO -(i/?i)(e*-e*-i 7 )(t-tt) „(0) 

~ ih 


99,P 


P^glS) = f dt'Epit'y 

J — oo 


Hh 

1 




“ pig,titp) 

= (4)2|n|2^0^0*g(i/fi)(e(-e‘+i7)(tp-t*)^W* 


99,P 


^(3),A(^) _ f^pF^Q*^-(i/h)(el-el-i-f)t^(i/h)eputp 
rt 


■ f dt'Ep{t')e^^/'^Fe-4-i'r-<^p-P'p^^J.{t') 

J —OO 


1 ^ ^0 (^i/n){el-el-ij)it-tp) (2) 






ih 

1 

^ih 

-(i/;i)(e*-e^-i7)(t-tp) (i/a)(e*-e*+i7)(tp-tt) (0)* 

^99,P' 
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Redefining the times as t — tp ^ t and tf — tp ^ r and 
recalling Eq. the signal contributing from the path A 
is given by 

Similarly, the signal associated with the path B and C 
are calculated as 

and 

In two-quantum regime (the positive delay), there are 
two pathways H and 1. We directly integrate Eq. [^and 
obtain the third-order polarization in the following way. 

P^egpi^) = 

J —oo 

in 

J —OO 

J —OO 

.p-b/^)(C-e?-*7)(tt-tp) AO) 

^ P99,P' 


With a redefinition the times, t — tf t and U — tp |r | 
and Eq. [^ the signal contributing from the path H is 
given by 

In the same way, the signal associated with the path I is 
calculated as 


where two different exciton modes are coupled. Note 
that we can easily apply our calculation to the polari- 
ton system in the main text just by rewriting the state 
|e) and \m) respectively to \L) {\U)) and \2L) {\2U)). 
In this case , the ground state energy is set to be zero 
and the following replacement holds: ef ^ ^l{u){^ 2 )i 

efn e2L(2C/)(^2), 4 ^ ^ ^2L{2U){ki), 

4 ^ ^L{u){—kFWM) and ^ ^ 2 L{ 2 U){-'I^fwm)- The 
coupling between polaritons and classical electric fields 
outside a cavity is represented by the quasi-mode con- 
pling \^L(U)- 


Appendix C: SIMULATION BASED ON 
EXCITON-PHOTON BASIS (LOCAL MODE 
SIMULATION) 

In this appendix, we discuss the connection between 
exciton-exciton interaction Hamiltonian Eq. and the 
polariton-polariton interaction Eq. Here, we assume 
that there is no motional (kinetic) degree of freedom by 
neglecting the energy momentum dispersion (the kinetic 
term) of exciton and photon. Eirstly, in the exciton- 
photon basis, the linear term Hun and exciton-exciton 
interaction term Hint are 

Hun = ex'^l'^x + eci’l'tpc + ^i'lpl'lpc + i>li>x) (27) 


Hint — ( 28 ) 

In the polariton basis, we can diagonalize the linear term 
^0 as 

Hlin = ^L^oi’Li’L + (29) 

The exact expression of the exciton-exciton interaction 
term Hint ia polariton basis is the following: 



(al) 

+ \9o\C\^'PuPu'^u'>Pu 

(a2) 

+ 2gQ\X\^\C\^'ip\tp\jtpLtpu 

(a3) 


(bl) 

- ga\C\‘^X*C^lti,]jti,u^u 

(b2) 

- go\X\'^XC*^\j^l^L^L 

(b3) 

- 5o|C'pXC*V>J^V>^V>;7V’L 

(b4) 

+ \gt,X*^Cmii,{hh 

(cl) 

+ ^9oC*‘^X‘^'lplj1pllpL'lpL- 

(c2) 


With the aid of the double-sided Eeynman diagrams, it is (al)-(a3) are respectively lower and upper polaritons self 
not difficult to extend the present discussion to the case and cross-interactions employed in the main text. The 
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terms (bl)-(b4) couple lower and upper-polaritons de¬ 
pending on the number of polaritons. The last terms (cl) 
and (c 2 ) annihilate two upper (lower) polaritons and cre¬ 
ate two lower (upper) polaritons. These terms are called 
the Darling-Dennison coupling terms [2|. Now we will 
evaluate the energy of the Hamiltonian H = Hun - 1 - Hint 
using a standard stationary perturbation theory of quan¬ 
tum mechanics [32] . The linear term Hun is a unpertur- 
bated Hamiltonian and the eigenstate of Hun is written 
as |n,m), where n and m are respectively the number 
of lower and upper-polaritons. Hint is the perturbative 
Hamiltonian. The energy = {n^m\H\n^m) is eval¬ 

uated as 


E — 


+ {n,m\Hint\n,m) 


E 


(n' ,m')jti(^n,m) 


E° - E°, , 

-^n,m 


+ ...(30) 


The first term, E^^^, is an unperturbed energy defined 
as E^ ^ = {n^m\Ho\n^m) = eL,on€u,o'm. The sec¬ 
ond and third terms respectively represent a first and 
second-order perturbation of the energy correction. A 
simple calculation shows that the self (al)-(a 2 ) and cross¬ 
interaction terms (a3) have diagonal elements and con¬ 
tribute to the first-order perturbation of energy. The 
other terms (bl)-(c 2 ) contribute only to the second-order 
of the perturbation. This means that the self and cross¬ 
interaction terms can be considered as dominant terms in 
the perturbative regime. The explicit form of the energy 
is calculated as 




+ \9o\X\'^n{n - 1) 

(a’l) 

+ ^9o\Cfm(m- 1) 

(a’ 2 ) 

+ 2go\X\'^\C\‘^n ■ m 

(a’3) 

+ ^|X|®X|C|V(n + l)m 

(b’l) 

+ ^\CfX\C\^{n + l)(m - ifm 

(b’ 2 ) 

-^^\X\^\C\\m + l){n-lfn 

(b’3) 

— ^ C ® Xp(m + l)rin?n 

(a’4) 

+ ^ X ^ C ^(n + 2 )(n + l)(m - l)m 

(c’l) 

- + 2 )(m + l)(n - l)n 

(c’ 2 ) 

+ ... 



The first term is a linear term of the unperturbed energy 
En^m- The next three terms (a’l)-(a’3) are the first-order 
perturbation terms resulting from the LP-LP (UP-UP) 
self and LP-UP cross-interaction terms. The remaining 


£ L{-k,) £l{ 0) £L{-ki] 


£iv(-/c,) £u{0)£u{-k,) 



-3 -2 

(c) 2-Quantum 

£'Li-k,} £l{ 0) £L{-k,) 


£u{k,) 
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£iv(-/c,) £u{0)£u{-k,) 



- 3 - 2-10123 
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FIG. 14. Simulation based on exciton-photon (local mode) 
basis Gross-Pitaevskii-equations. Amplitude of delaymap 
|5'(t, et)| (a). Amplitude of 2D FWM spectrum |5'(er,et)| 
for 1-quantum (b) and 2-quantum region (c). 


parts are the second-order perturbation terms. In the 
main text, we neglect the second-order terms (b’l)-(c’ 2 ). 
Let us evaluate the condition where this approximation 
holds assuming that the lower and upper-polariton den¬ 
sities (m = n) are equal and the cavity detuning is 
zero (|X| = \C\). The calculated energy indicates that 
the first-order perturbation term is proportional to gon^ ^ 
while the second-order one is to gon^/Q. Thus, the 
approximation which takes into account the first three 
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terms is satisfied under the condition: /Vt < 1. 

Hence, a large Rabi splitting and low polariton density 
is required for our approximation to hold. Actually, if 
the number of polaritons becomes large, this perturba¬ 
tion breaks down and the contribution from the terms 
(bl)-(c2) becomes comparable to (al)-(a2). For sake of 
rigor, we note that our treatment only holds when the 
polariton-polariton scattering terms away from q=0 can 
be ignored. 

In addition to the above perturbative discussion, we 
present 2DFT spectra based on exciton-photon Hamil¬ 
tonian Eq. using numerical simulations. Similarly to 
the polariton Gross-Pitaevskii equations, with a mean- 
field approximation, we can derive conventional exciton- 
photon Gross-Pitaevskii equations [18]: 

ihiix = {(^x - + go\'ii’x? - i^)^x (31) 

ih-lpc = (Cc - -i^)'>Pc- /ext- (32) 

Irric 2 


The simulation of the cavity detuning ^=-0.38 meV is 
presented in Fig. The parameters used in this sim¬ 
ulation is same as that of polariton-based calculation 
performed in the main text (Fig. [^. In this exciton- 
photon basis calculation, both nonparabolicity of po- 
lariton’s energy-momentum dispersion and all polariton- 
polariton interaction terms Eq (al)-(c2) are automat¬ 
ically taken into account. The comparison between 
exciton-photon (Eig. 14) and polariton basis calcula¬ 
tions (Eig. show that both frameworks give quali¬ 
tatively the same results for low polariton densities. In 
the high density regime, the results of the two simula¬ 
tions change (not shown), where the lower-polariton basis 
calculation breaks down, while the exciton-photon basis 
Gross-Pitaevskii still works. 
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